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Motivated by the nickel valence controversy in the perovskite nickelate RNi03, we consider a 
one-dimensional underscreened Kondo chain consisting of alternating spin-1 ("nickel") and electron 
("oxygen") sites, which in addition to the usual electron hopping and spin-spin interaction between 
the 5 = 1 spin and the electron also contains a spin mediated electron hopping term. Using the 
density-matrix renormalization group (DMRG), we obtained the zero temperature phase diagram 
of the model, as well as various correlation functions in each phase. Importantly, for certain range 
of parameters the model exhibits a quasi-long-range spiral (QS) order. To understand the DMRG 
results, we construct a mean-field theory based on a Schwinger fermion decomposition of the S = 1 
spins, from which we argue that the QS phase corresponds to a phase in proximity to the spin Bose 
metal state proposed by Sheng, Motrunich, and Fisher [Phys. Rev. B, 79, 205112 (2009)]. Notably, 
we find no evidence for a phase with the symmetry of "nickef'-centered charge order, which has 
been argued to arise due to site-selective Kondo screening of half the S = 1 spins, and suggest that 
order of this type occurs only due to an additional energy gain from spontaneous lattice distortions, 
not present in this model. 



I. INTRODUCTION 

Mott metal-insulator transitions and their associated 
spin and charge ordering appear throughout the physics 
of transition metal oxides, and present many puzzles 
of interpretation and modeling. Since electrons in the 
regime of a Mott transition may possess both localized 
and itinerant character, ambiguity may arise as to an ap- 
propriate model Hamiltonian: Hubbard, Heisenberg, or 
t-Jl Anderson or Kondo model? In multi-orbital sys- 
tems, the choices are even more numerous. Such a sit- 
uation arises in the pseudo-cubic perovskite nickelates, 
RNiC>3, where R is a rare earth ion, which have been 
studied for many years for their interesting Mott metal- 
insulator transitions. 1 - A T > transition occurs to an 
insulating state for all members of the series save the 
case R=La, which is metallic at all temperatures, to a 
low temperature state with complex structural symme- 
try breaking — denoted "charge order" in the literature — 
and antiferromagnetism. Nominally in this material the 
nickel valence is 3+, i.e. 3d 7 , which would suggest a 
model with a single electron in the cubic e g doublet, a 
two-fold orbital degeneracy. However, another interpre- 
tation is that the nickel valence is actually 2+, with an 
additional hole per nickel spread out amongst the oxygen 
ions, analogous to the Zhang-Rice singlet in cuprates^. In 
this "ligand hole" scenario there is no orbital degeneracy, 
but rather an S = 1 spin on the nickel site, correspond- 
ing to a half-filled e g doublet. The difference between 
these two scenarios is well-defined only in the ionic limit 
where the charge on the nickel site does not fluctuate 
significantly, which might be the case in the ligand-hole 
picture. A third, compromise view is that the valence 
fluctuates significantly due to hybridization, making the 
distinction between Ni 2+ and Ni 3+ moot; in this case one 
might build a theory in terms of the low energy bands 



near the Fermi energy. Though such a compromise exists, 
it seems to be unpopular, and strong opinions in favor 
of the ligand-hole picture are often voiced. Going fur- 
ther, Sawatzky has proposed a picture for the observed 
"charge" ordering in terms of local Kondo screening of 
half the Ni 2+ spins. 

In this paper, in hopes of gauging the plausibility 
of the ionic ligand-hole scenario, we study a concrete 
model using numerically exact methods: the Density Ma- 
trix Renormalization Group (DMRG). 3 ' We write down 
the simplest model incorporating a sharply-defined Ni 3+ 
valence, and holes on the oxygen sites. Because each 
nickel thereby has an S = 1 spin and is accompanied 
by only a single oxygen hole with S = 1/2, this takes 
the form of an underscreened Kondo lattice modelW^ To 
enable the DMRG analysis, we take this model to be 
one-dimensional. Though this is obviously a drastic ap- 
proximation, our model retains the underscreened Kondo 
physics, and has the symmetries to allow charge ordering 
of the same type observed in the nickelate materials. We 
obtain the numerically exact full zero temperature phase 
diagram for this model, which contains several magnetic 
states. However, we do not find the charge ordered state 
seen in the nickelates, which is nickel site-centered, but 
instead a complementary type of charge ordering which is 
bond- centered, i.e. ordered on the oxygen sites. We will 
discuss the implications of this result for the nickelates 
in Sec. EH 

Beyond the context of the nickelates, our model is 
of interest on its own as a problem in one-dimensional 
physics, and of the underscreened Kondo effect. Recent 
work has shown the possibility of gapless one-dimensional 
"spin liquid" phases, which are strikingly even more non- 
quasiparticle in nature than the usual Luttinger liquid. 
Such phases, have been found in rather exotic models 
with very large ring exchange interactions.^ Remarkably, 
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FIG. 1: (color online) The nickelate chain and the t-J-J' 
Hamiltonian (bottom panel), as derived from a Hubbard- 
typed model (top panel). 



we find evidence for such a spin liquid-like phase in our 
rather simple and physically motivated Kondo lattice 
model (technically, we obtain a phase which appears to 
be proximate to the "spin Bose metal" phase of Refill 
rather than this phase itself) . This suggests an intriguing 
prospect of observing quantum spin liquid states in two 
and three dimensional Kondo lattice systems appropriate 
to real materials. Again, this idea has already been sug- 
gested in the literature^, but connections do not seem 
to have been made to realistic model Hamiltonians. 

The remainder of this paper is organized as follows. In 
Sec. |TI] we introduce our model Hamiltonian as well as 
a simple semi-classical treatment that would guide our 
intuition about the quantum case. In Sec. |III| we present 
our DMRG results, which includes the quantum ground- 
state phase diagram and various correlation function in 
each phase. In Sec. |IV| we introduce a mean-field pic- 
ture, which, together with degenerate perturbation the- 
ory, explains most phases in the DMRG phase diagram. 
In Sec. [V] we focus on the one phase that cannot be ex- 
plained in the previous section, which we argue is a phase 
in proximity to the spin Bose metal phase in Ref. H3 Fur- 
ther discussions and conclusions are presented in Sec.|VI| 



II. THE t-J-J' HAMILTONIAN 

Motivated by the nickelate, we consider a one- 
dimensional (Id) chain consisting of alternating electron 
and spin-1 sites, as illustrated in Fig. [TJ in which the 
electrons are half-filled. We shall continue to refer to 
the electron sites as the oxygen sites and the spin-1 sites 
as the nickel sites, even though our model is no longer 
constrained by material details. 

The interactions between the electrons and the nickel 
spins arise from virtual hoppings. At second order, after 
hopping to the nickel site, the electron or hole can either 
hop back to the original oxygen site or to the next oxy- 
gen site, which gives rise to two distinct contributions to 
the Hamiltonian. In addition, the electron can also hop 
directly from one oxygen site to the next without the 
mediation of the nickel spin in-between. Taken together, 



this leads to the t-J-J' Hamiltonian: 
Ht J J' — — t c 1 . ! c. ! + h.c. 
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where c. , , and , , are the electron operators on the 

oxygen site at i ± \ with spin index a, Sj is the S = 1 
spin operator on the nickel site i, and <x is the vector of 
the usual Pauli sigma matrices. H t , Hj, and Hj* are 
defined in the obvious way. For convenience, we also 
define s,- + i = c' , ^W-c, i „. i.e., s,- + i is the electron 

spin operator on site i±s- Note that here and henceforth 
the spin indices are assumed to be appropriately summed. 

It is worth noting that the somewhat unfamiliar J' 
term in H t jji, is analogous to the density dependent 
hopping term that is customarily neglected in the t-J 
model. 1 -- However, unlike in the t-J model, for which at 
low-doping this term can be approximated by a density 
independent hopping, in the present case the spin vari- 
able on the nickel site remains strongly fluctuating even 
at low energy and thus the J' term cannot be neglected. 
This will become evident in Sec. Mil 

It should also be noted that J and J' are in general re- 
lated to each other through the parameters of the under- 
lying Hubbard-type Hamiltonian (in which charge fluc- 
tuation is put back to the nickel sites) and in the most 
natural scenario J = J' . However, in this manuscript 
we will treat J and J' as individual parameters without 
worrying about how they may arise from a Hubbard-type 
Hamiltonian. We shall, however, restrict ourselves to the 
quadrant in which J, J' > 0. 

To gain some intuitions about the t-J-J' model, con- 
sider a semi-classical treatment in which the nickel spins 
are treated as classical while the electrons remain quan- 
tum mechanical. For any configuration of nickel spins, 
Htjj' reduces to a quadratic Hamiltonian of the elec- 
trons at half filling, and the ground-state configuration 
of the nickel spins is determined by minimizing (H t jj>) 
with respect to all possible nickel spin configurations. 

For simplicity we consider only spiral configurations 
of the form Sj = S (cos(qXi)x + sin(qxi)y) for the 
nickel spins. In such case, the diagonalization of the 
quadratic Hamiltonian is facilitated by the transforma- 
tion c. . i — > e ±lqXi l 2 c. , i , where the + (— ) sign holds 

for a =f (4-). Such transformation removes the position 
dependence in the coefficients of the quadratic Hamilto- 
nian obtained from H t jj', which in turns can be diago- 
nalized by a simple Fourier transform. 

Carrying out the minimization with respect to q, we 
find the phase diagram shown in Fig. [2j in which the 
nickel spins form an antiferromagnet (AF; q — ir) when 
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FIG. 2: (color online) Phase diagram obtained from the 
semi-classical computation, where AF stands for a antiferro- 
magnetic state in which wavevector q = n, F stands for a 
ferromagnetic phase in which wavevector q — 0, and S stands 
for a spiral phase in which < q < n. Note that all three 
phases meet at the point J — J' — t. 



J' 3> J, a ferromagnet (F; q = 0) when J ^> J', and a 
spiral (S; < q < n) state when J m J' > t. 

However, it is well-known that under general circum- 
stances the continuous symmetries of a Id quantum sys- 
tem cannot be spontaneously broken. In particular, in 
the well-known antiferromagnetic J\-Ji model, the spi- 
ral phase of the classical model turns into a dimer phase 
when the classical spins are replaced by quantum S = 1/2 
spins. In the dimer phase, the spin-spin correlation 
(s_fc • Sfe) exhibits no singularities but only peaks, first 
at q — 7r and then at an incommensurate wavevector 
as J2 increases! 11 * 12 ^ With this in mind, it is natural to 
ponder about the fate of the spiral phase in the present 
model when the quantum mechanical effects are taken 
into account. 



III. DMRG COMPUTATIONS AND RESULTS 

We determine the quantum ground-state phase dia- 
gram of the t-J-J' Hamiltonian, Eq. ([I]), by large-scale 
density-matrix renormalization group (DMRG^P calcula- 
tions. We consider systems with up to N = 96 unit cells, 
which consists of oxygen and one nickel site (see Fig. [3] 
for illustration). Specially, in the DMRG calculation, we 
use open boundary conditions (OBC), and keep up to 
to = 4000 states in each DMRG block. This is found 
to give excellent convergence in the measurements such 
as the ground state energy and various correlation func- 
tions with total error of the order or less than 10 -6 . The 
phase boundaries in the {J/t, J' /t) parameter space are 
determined by extensive scans of the derivatives of the 
ground state energy and by monitoring the correlation 
functions, as well as the corresponding order parameters. 

To determine the properties of the ground states, 
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FIG. 3: (color online) DMRG block used in the calculations. 



TABLE I: Charge gap A and central charge c for characteris- 
tic parameter values within each phase of the phase diagram 
Fig. [4] For the CD phase, "com." indicates peak in (S k Sl k ) 
at k = 7r while "incom." indicates peak in (S k Sl k ) at incom- 
mensurate wavevectors. 
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we calculate the nickel spin-spin correlation (5'^S'° fc ), 
the electron spin-spin correlation (s^s^), the elec- 
tron density (n^), the electron density-density corre- 
lation (dnkSri-k) (Sn is defined in position space by 
5rii = tii — 1), and the nearest-neighbor nickel-oxygen 

spin-spin correlation ^s i± i ■ Sj^. To determine if the 

phase is conducting, we also calculate the charge gap 
A = [E Q (N + 2) + E (N - 2) - 2E (N)}/2. Finally, we 
also obtain the central charge c by calculating the von 
Neumann entanglement entropy Sa for a partition of the 
system into two halves of length i, L — t, with varying 
length I of the sub-system A. A universal logarithmic de- 
pendence of Sa on £ is expected in the thermodynamic 
limit with a prefactor depending upon c, and accurate re- 
sults can usually be obtained using finite-size scaling for- 
mulae derived from conformal field theory, see Refs ll3ll4"l 
The calculation of the von Neumann entropy is compu- 
tational expensive and generally the most difficult of all 
physical quantities to converge, especially in highly en- 
tangled states. Therefore our results for c are more lim- 
ited than for the other quantities discussed above. 

Our main result is the phase diagram presented in 
Fig. [4| in which we find four distinct phases as J and J' 
are varied. In addition, the correlation functions at char- 
acteristic points in parameter space are shown in Fig. [5] 
and [6j Note that we have only shown the zz-components 
of the spin-spin correlation functions. We remark that 
the phases labeled as QAF, CD, and QS (detailed below) 
possess spin SU(2) invariance, and thus the other com- 
ponents of the spin-spin correlations are essentially iden- 
tical to the ones we showed. For the ferromagnetic (F) 
phase, the z-axis in spin space is constrained to M z = 0, 
and thus the magnetization axis is perpendicular to the 
z-axis. 

For J > J' the system exhibits ferromagnetic order 
with magnetization M ~ 0.5 (1b P er unit cell. This 
is consistent with previous works in the literature on 




FIG. 4: (color online) Phase diagram from the DMRG cal- 
culation. Here F stands for the ferromagnetic phase, QAF 
stands for the quasi-long-range antiferromagnetic phase, CD 
stands for the charge-density ordered phase, and QS stands 
for the quasi-long-range spiral phase. The dashed line within 
the CD phase represent boundary in which the peak in 
{SlS*L k } changes from fc = 7rtofc = 7r±<5. 



the alternating spin- 1-spin- 1/2 chainp^E^D which can be 
thought of as the J' = 0, J — > oo limit of the present 
model. For J and J' both small and J' > J, the system 
exhibits quasi-long-range antiferromagnetic (QAF) order 
characterized by sharp peak in (S k S^_ k ) at k — n, while 
( s k s -k) an d (SntSn-k) remain largely featureless. As J' 
further increases such that J' > J w t, the system devel- 
ops a charge-density (CD) order with period 2 as shown 
in the inset (a) in Fig. [7j and the sharp singularity in 
(S k S^_ k ) in the QAF phase becomes a broad peak. As 
J increases within this CD phase, the peak originally lo- 
cated at k = n splits into two peaks at k — tt ± S, with 5 
increasing as J increases. 

Most interestingly, when J and J' are both large and 
J' > J the system seems to exhibit quasi-long-range spi- 
ral (QS) order characterized by sharp singularities in both 
the nickel spin-spin correlation (S k S*L k ) and the electron 
spin-spin correlations (s^s z _ k /j at incommensurate vec- 
tors k — 7T ± 5, while the electron density correlation 
(SnkSn-k) remains largely featureless, and no signs of 
CD order are found. Moreover, when the correlation 
functions are examined more carefully near k = n, a 
tiny broad peak can be discerned at wavevector roughly 
equal to ±<5 in (S k S^_ k ), while no such features are found 
in (s£sf_ fe \. It should be remarked that the value of 5 
evolves continuously as J' and J changes and connects 
across the phase boundaries into the ferromagnetic and 
the CD phases, as shown in Fig. [8] 

In addition, from Fig. [6] we see that in the QAF and 
the QS phases, the nearest-neighbor spin-spin correlation 

^s i± i • Sj\ is essentially uniform in the thermodynamic 

limit if not already in finite-size systems, while in the 
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FIG. 5: (color online) (a) Nickel spin-spin correlation 
{SkS-k)i 00 electron spin-spin correlation {s%s*L k }, and (c) 
electron density-density correlation {5nkSn-k) for (orange, 
filled circular symbols) J/t = 0.4 and J'/t = 1.6, (green, 
empty square symbols) J/t = 0.8 and J'/t = 4.0, (blue, 
empty diamond symbols) J/t — 1.6 and J'/t = 4.0, and (pur- 
ple, greek cross symbols) J/t — 2.4 and J'/t — 4.0. The 
insets of panels (a) and (b) show the details of the respec- 
tive correlation for J/t = 2.4 and J'/t = 4.0 in the range 
0.5 < k/TY < 1.5. 



CD phase \s i± i ■ Sij has a significant two- unit-cell os- 
cillation centered at the oxygen site, which is consistent 
with the symmetry of the CD order. Note that the os- 
cillations of ^Sj_|_i ■ Si) in the CD phase have larger am- 
plitudes than the corresponding oscillations in (n,) , and 



that in general 



■Si 



increases as the parameters 



(J/t, J'/t) increases. 

To shed further light into the phases we obtained in 
DMRG, we show the charge gap A and central charge c 
for characteristic parameter values within each phase in 
Tab le |T| (computational details are illustrated in Figs. [9] 
and 10). From the table it can be seen that all phases 
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FIG. 6: (color online) The nearest-neighbor nickel-oxygen 

spin-spin correlation (s i± i ■ Sjy for the characteristic points 

in parameter space used in Table [I] with system size N — 48. 
The two correlations are combined and the index i now labels 
the bond center of the correlation. The positions of the nickel 
and the oxygen sites are indicated by the labels above the 
curves. The inset shows finite size scaling of the oscillation 

amplitude in /s i± i • Si) for J = 0.4 and J' = 1.6 using a 

second order polynomial fit. 



have finite charge gaps. Moreover, the values of the cen- 
tral charges suggest that one bosonic degree of freedom 
remains gapless in the QAF phase, while no degree of 
freedom remains gapless in the CD and F phase. Unfor- 
tunately, owing to the large (twelve) internal dimension 
of the unit cell, we cannot accurately determine the cen- 
tral charge of the QS phase. 



IV. MEAN-FIELD KONDO PICTURE AND 
EFFECTIVE SPIN-1/2 MODEL 
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FIG. 7: (color online) The development of charge-density 
wave order parameter Sncu = {tin) — (njv +1 ) as J'/t in- 
creases while J/t = 0.8 is fixed. The inset (a) shows (m) as 
function of unit-cell index i for (black, filled square symbols) 
J/t — 0.8 and J'/t — 2.4 and (blue, filled diamond symbols) 
J/t = 0.8 and J'/t = 4.0, with system size N = 96. The Inset 
(b) shows examples of finite-size scaling of <5ticd at different 
J' using second order polynomial function. 




FIG. 8: (color online) Locations of the singularity (resp. 
peak) in the ferromagnetic and QS phases (resp. CD phase) 
as function of J'/t when J/t — 2.0 is fixed. 



In this section we introduce a slave-particle representa- 
tion of the nickel spin, from which we obtain a quadratic 
mean-field theory (subsection IV A). Assuming the sim- 
plest form of mean-field ansatz, this leads to a spectrum 
that contains a flat band exactly at the Fermi energy. 
We then show that such a degeneracy can be lifted us- 
ing second-order degenerate perturbation theory, which 
allows us to map the present problem to an effective spin- 
1 /2 model ( IV B I . We then carry out an explicit calcu- 
lation of the parameters in the effective spin-1/2 model, 
and compare a quantum mean-field analysis of the phase 
diag ram o f this model with the DMRG results (subsec- 
tion IV C I . Our main result in this section is the quantum 



mean-field phase diagram of Fig. 16 



A. Slave-particle representation and its mean-field 

To interpret the DMRG results, it is useful to intro- 
duce a mean-field picture. Specifically, we want mean- 
field "parent states" that are translationally and spin 
SU(2) invariant, under which the quasi-long-range orders 
and/or the broken symmetries manifest as singularities 
in the correlation functions in the resulting low-energy 
effective theory. To this end, we consider the following 
Schwinger fermion representation, in which the spin-1 op- 
erators Si are written in terms of four species of spinons 
f iaa , carrying both orbital indices (indicated by lowercase 
Latin characters) and spin indices (indicated by lowercase 
Greek characters): 
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FIG. 9: (color online) Finite size scaling of the charge gap for 
the characteristic points in parameter space used in Table [I] 
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FIG. 10: (color online) The lower branch of the von Neu- 
mann entanglement entropies, obtained from the lower enve- 
lope of the full entropy shown in the inset, as a function of the 
conformally-transformed length i' of the subsystem A for the 
characteristic points in parameter space in Table [I] The solid 
lines are linear fits using Sa = § ln(i'), which determines the 
central charge c. 



Substituting the Schwinger fermion representation into 
the t-J-J' Hamiltonian results in four-fermion terms that 
are schematically of the form c^f'fc, which can be han- 
dled by a Hartree-Fock-type decomposition. Together, 
this leads us to the following mean-field Hamiltonian: 



ij MF = - Me w c -<y c f +1 c. , +h. c . 

i i 

i 

+E A /L/ ia a + N-/L^ ab / lbQ , (4) 



in which the Lagrange multipliers A and N are introduced 
to enforce the constraints Eq. ([3| on average, and that the 
orbital-dependent Kondo hoppings Q a and Q' a are to be 
determined by the minimization of the mean-field energy, 
i.e., Q a and Q' a are to be chosen such that (H t jj<) MF 
is minimized (here (-) MF denotes expectation value with 
respect to the mean-field state obtained from H MF ). 

Notice that A, N, Q a , and Q' a carry no site indices 
since translation invariance is assumed. Similarly, Q a 
and Q' a carry no spin indices since spin SU{2) invari- 
ance is assumed. By the same token, we also ne- 



glect the "magnetic" mean-field terms {f} aa & ati f ia p) and 

(ct, x <T afi c. , J mH MF . 

It should be remarked that H MF can alternatively be 
derived in a Feynman path integral approach in which 
fluctuating A^ and Nj are introduced as auxiliary fields 
to the partition function, such that the constraints are 
enforced exactly upon functional integration. Similarly, 
in this approach Qi a and Q' ia are introduced as fluc- 
tuating bosonic Hubbard-Stratonovich fields that upon 
functional integration reproduces the appropriate four- 
fermion termsPSMl i n this context, Ai and can also 
be interpreted as the temporal components of a U(l) and 
an SU (2) gauge field, respectively. Together, they corre- 
sponds to the U(2) ~ U(l) x SU(2) gauge redundancy 
fiaa U°" b fib a of the spinons. However, at this stage we 
shall take the simpler picture and treat the variables A, 
N, Q a , and Q' a as parameters of H MF . 



a/3 
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in which the spinons are subjected to the constraints 
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E/L^ Qb / 4hQ =0, (3) 
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(recall that the spin indices are implicitly summed), in 
which t is the vector of the Pauli sigma matrices that 
act on the orbital indices.^ For brevity, henceforth we 
shall assume that the orbital indices are appropriately 
summed. 



B. Mean-field ansatze, flat bands, and degenerate 
perturbation theory 

While in general Q a and Q' a are distinct, it is natu- 
ral to first consider scenarios in which Q a = Q' a , where 
we can make use of the U(2) gauge redundancy to set 
Qa = Qa = r[l,0] T , with r > 0. Then, the con- 
straints Eq. ^ requires N oc z and thus the last line 
of Eq. fl can be rewritten as J2 ia ^a,f} aol f iaa , in which 

\ a = \+(-l) a+1 N z . 

From this, it can be seen that for this class of ansatze, 
the ji spinons are completely decoupled from the other 
fermions in the system and enter H MF only through a 
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FIG. 11: Dependence of r/t on (J + J')/t in the class of 
mean-field ansatze denned by Q a = Q' a = r[l,0] T . 
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FIG. 12: Typical mean-field spectra obtained from // for 
the class of mean- field ansatze defined by Q a = Q' a = r[l, 0] T . 
Note that the degenerate flat bands at zero energy are not 
plotted. 



(specics-dcpcndcnt) chemical potential. Consequently, 
the f% spinons form a flat band in the mean-field spec- 
trum, which has to be half-filled in order to satisfy the 
constraints. The two remaining species of fermions form 
a Kondo band insulator, in which the band gap is con- 
trolled by the ratio r/t that increases as J, J' increases. 
Moreover, it is easily checked that for this class of ansatze 
(Hj) = (Hj>), and hence the mean-field spectrum de- 
pends only on (J + J')/t. Carrying out the minimiza- 
tion of (H t jjr) MF with respect to r, we obtain r/t as 
a function of (J + J')/t as shown in Fig. [TTJ The 
mean- field spectra for ansatze with different r/t are plot- 
ted in Fig. [12] Note that in the mean-field picture 



oc 



t ax: 



Hence, the 



trend of increasing r/t as (J + J')/t increases is consis- 
tent with the trend in DMRG [c.f. Fig. [6]. 

Physically, the half-filled flat band can be interpreted 
as free 5 = 1/2 spins on the nickel sites that arise from 
Kondo under-screening, in which the electrons screen out 
only half of the nickel spin on each site (see Fig. 13 for 
illustration). This degeneracy among the ]i spinons is 



FIG. 13: (color online) Schematic of the Kondo under- 
screening, in which the electrons screen out only half of the 
nickel spin on each nickel site, leaving behind an unscreened 
spin- 1/2. 



expected to be lifted when effects beyond mean-field are 
considered. To capture such effects, we perform second- 
order degenerate perturbation theory on the mean-field 
state, in which the perturbation is provided by the resid- 
ual interaction H res = JHj + J'Hj, - r ^Ac] 1 Jiia + 

1 f ;l „ + h.c). Notice that the constraint terms do 
not enter H les , since the constraints will still be satisfied 
after the degeneracy of the fi spinons are lifted. 

As usual, from second-order degenerate perturbation 
theory we obtain an effective Hamiltonian H , given by: 



H 



err 



V G H res V x 



1 



£ - H MF 



V x H Ies V G , 



(5) 



where Vg is the projection operator onto the degener- 
ate ground-state manifold, Vx = 1 — Vg, and £q is the 
unperturbed mean-field ground-state energy. 

Since the f% spinons are decoupled from the electrons 
and the /1 spinons, the eigenstates \ip) of H MF can be 
written as product states of the form \%b) c f 1 &> \ip)f 2 , an d 
correspondingly the Hilbert space decomposes as J4? = 
J4? c f! < £>J4°f 2 - In this language, the ground-state manifold 
of H MF is {IG)^} (g) Jij 2 , where |G) C / 1 is a unique non- 
degenerate state. Moreover, it can be checked that H res 
factors into the following form: 



H - I VT ■ f f ' 



■a8 



-I 



2i8 



(6) 



where Tj and the ". . ." operates within -^cfx- Explicitly, 



t I f & f n 

" I C 1 ~ C . , 1 n — 1 C . i T C . In 



+ J' ( CI 1 — — C . 1 a + C 

Using Eq. (joj), H cS becomes: 
7- I f f - 

J n I Jila 



r Q/3 



r Q/3 



(7) 



= JijSi ■ Sj + const. , 




const. 



(8) 



where from translation symmetry it follows that Jij de- 
pends only on \xi — Xj\, i.e., Jij — J\i~j\- Recognizing 

fi2a^~T~ fi2/3 = s i as a spin- 1/2, the second-order degen- 
erate perturbation theory thus maps the t-J-J' model to 
an effective spin- 1/2 model. 



Classically, the ground state of the effective spin model 
Eq. d8) is a spiral state given by Sj = cos(xj9)n-i + 
sin(xJ9)n2, where fii and 112 are two orthogonal unit 
vectors, and which 9 is chosen to minimize the classical 
energy. However, in Id, this classical order is expected 
to be destroyed by quantum fluctuations. 



C. Explicit calculation of Jij and comparison with 
DMRG results 

Returning to the present case, J7ij in Eq. ([8| equation 
is given schematically by: 



J h 



c h (G\T?\X) ch cfl (X\T*\G) cfl 



+ c.c. 



(9) 

in which \X) c f 1 are particle-hole excitations from |G) C / 1 , 
with £x its energy. Note that we have made use of the 
spin 517(2) symmetry to evaluate Jij using only the z- 
component of TV 

To write down Jij more explicitly, we Fourier- 
transform and diagonalize the part of H MF that involves 
only the electron and the /i spinon, H^jF, as follows: 



rrMF 



EE 

k At=± 



(10) 



where by convention e + k > e_^. We also define the 
eigenvectors ui via the following: 



(11) 
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Cka 




'< + 






fka 




/+ /- 




l-,k,a 



Plugging in, we arrive at: 

1 >r-^ cos ((k — q)(Xj — Xi)) 



Jij — 



k.q 



e -,k - ^ 



x \ j{l + e Lq e- lk ) + J'{e lq + e- lk )\ 



(12) 



where N is the number of unit cells. 

Importantly, the J and J' terms in the above equation 
carry different Fourier phase factors. Thus, even though 
the ratio r/t in the mean-field ansatze depends only on 
(,/ + J')/t, different effective spin-1/2 models can still be 
realized for different J and J' having the same sum. 

Calculating Jij from Eq. ( 12 ) up to the eighth near- 
est neighbor, we obtain the classical phase diagrams as 
shown in Fig. |14| in which we distinguish between the an- 
tiferromagnetic phase (9 = tt), the ferromagnetic phase 
(9 = 0), and the spiral phase (0 < 9 < it). We also plot 
the classical spiral angle 9 as function of J' /t along the 
line cut J' — 2 J and J' = 4 J in Fig. 
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From these figures it can be seen that even though the 
details of the phase boundary and the spiral angles are 
modified as further-neighbor interactions are included, 
the truncation at Ji still captures the qualitative as- 
pects of the model reasonably well. Hence, we now focus 
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FIG. 14: (color online) Classical phase diagram of the ef- 
fective spin-1/2 model Eq. (|8| truncated at (a) Ji and (b) 
Js- Here filled gray triangular symbols indicate the ferro- 
magnetic phase, filled orange circular symbols indicate the 
antiferromagnetic phase, and unfilled teal square symbols in- 
dicate the incommensurate spiral phase. The lower-left corner 
of the phase diagrams are excluded because the gap in H MF 
is too small to accurately calculate Jij. 



on the results obtained within this truncation and map 
the parameters J\ and J2 we obtained to the known re- 
sults from the quantum J\-Ji model, in which the ground 
state is known to exhibit QAF order when Ji ^S> | J2I > 0, 
and undergoes a quantum phase transition into a dimcr 
state at J2/J1 ~ 0.241. Moreover, as J2 further increases 
beyond J2/J1 = 1/2, the peak in the spin-spin correla- 
tion originally located at k = 7r also splits into two in- 
commensurate peaks at k = 7r ± S, with S increasing as 
d = tan _1 ( J2/J1) increases. Eventually, the system b e- 
comes ferromagnetic when ■& > tan -1 (— l/4).mi 2 l 21 l 22 l. 

The results of our mapping from the (J/t, J/t) pa- 
rameter space to the J1-J2 model is shown in Fig. |16| 
Note that in the present case, there are two inequivalent 
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FIG. 17: (color online) Illustration of the two inequivalent 
oxygen sites when the effective spin- 1/2 model is in the dimer 
phase. 



FIG. 15: (color online) Classical spiral angle 6 of the effective 
spin- 1/2 model Eq. Q as function of J' along the line cut (a) 
J' = 4J and (b) J' = 2J. The red, solid (blue, dashed) curve 
with filled circle (cross) symbols corresponds to truncation at 
J2 (Ja). 
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FIG. 16: (color online) Quantum phase diagram of the effec- 
tive spin-1/2 model Eq. jjjj, with Jij computed from Eq. ( 12 1 
and truncated at J%. Here filled gray triangular symbols in- 
dicate the ferromagnetic phase, filled orange circular sym- 
bols indicate the quasi-long-range antiferromagnetic phase, 
unfilled green square symbols indicate the dimer phase with 
spin-spin correlation peaked at k — ir, and unfilled violet dia- 
mond symbols indicate the dimer phase with spin-spin corre- 
lation peaked at an incommensurate wavevectors. The lower- 
left corner of the phase diagram is excluded because the gap 
in H MF is too small to accurately calculate Jij. 



oxygen sites when the effective spin-1/2 model is in the 
dimer phase, as illustrated in Fig. [17] From this, it can 
be seen that the broken symmetry in the dimer phase of 
the effective spin-1/2 model is precisely the broken sym- 
metry one would expect from the period-2 CD phase of the 
DMRG phase diagram (Fig. [^J). Moreover, even though 
there is no charge degree of freedom left in the effective 



spin-1/2 model, charge order is likely to occur when fur- 
ther effects beyond mean-field, e.g. the back-reaction of 
the dimer order onto the Kondo band insulator formed 
by the electron and the f\ spinon, are taken into account. 
Given that the charge deviation in the DMRG CD phase 
is small, this scenario in which charge order is derived 
from the ordering of the spin degree of freedom is con- 
sistent with the DMRG results. Therefore, we identify 
the dimer phase of the effective spin-1/2 model with the 
DMRG CD phase. 

Note that in this picture, the charge order in the CD 
phase is driven by the spin order. Thus, we expect a 
larger oscillations in spin correlations than in charge den- 
sity, consistent with the DMRG results [c.f. Figs. ^ and 
nh]. We also remark that the amplitude oscillation in 

(^ c l±i/2-fii) ~ \]\( s i±i/2 ' Si) | cannot be obtained at 
the mean-field level even if we extend the unit cell to 
two nickel and two oxygen per cell and allow r to vary 
from bond to bond, as long as the degeneracy of the fi 
spinons are left untouched. This highlights the impor- 
tance of the dimer formation in the effective spin-1/2 
model as the driving mechanism of the oxygen-centered 
dimer/CD order. 

Combining, we see that Fig. [16] captures the essential 
aspects of the DMRG phase diagram well, except for the 
QS phase. This is particularly so if one allows for sepa- 
rate renormalizations of J and J' from their bare values, 
which one can easily imagine to have occurred when vari- 
ous effects that we have neglected are taken into account. 



V. MEAN-FIELD IDENTIFICATION AND 
EFFECTIVE THEORY OF THE QS PHASE 

In this section we focus on the QS phase by extending 
the class of the mean-field ansatze considered. Some- 
what surprisingly, the flat bands persist even after the 
class of ansatze under considerations is extended (sub- 
section 



V A I . To lift such degeneracies without excessive 



complications, we introduce an additional slave-fermion 
hopping parameter t' , from which we obtain modified 
mean-field spectra in which for generic parameters a sin- 
gle band crosses the Fermi energy at two pairs of Fermi 
points, (subsection VB). Next we bosonize the effective 
theory obtained from the low energy fermions (subsection 
V C ) and search for an appropriate combination of effec- 



tive interactions that best reproduces the salient features 
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FIG. 18: (color online) Results from energy mini- 

mization for the general mean-field ansatze. Here un- 
filled blue circle (red cross) indicates states for which 

\Jl2a\{fL^ a ) ~ (/L^_1 Q )| 2 < (» 0.05. 



of the QS phase (subsection V D I 



Our main result in this section is the identification of 
the QS phase with a bosonized theory that has one spin 
and one charge ("C1S1") mode, in which there is a finite 
charge gap and the spin field carries an incommensurate 
wavevector. This interacting bosonized theory is prox- 
imate to the spin Bose metal state proposed by Sheng, 
Motrunich, and FisherP 



A. Extended mean-field ansatze 

Since the class mean-field ansatze restricted to Q a = 
Q' a fail to describe the QS phase, we now consider general 
ansatze in which Q a and Q' a may be unequal. As in the 
preceding section, the U{2) gauge redundancy f iaa — > 
U^ b fiba can be used to reduce the number of parameters 
needed to specify all physically distinct ansatze. To be- 
gin with, it is clear that we can fix Q a = r[l,0] T , with 
r > 0. However, such choice does not exhaust the U(2) 
gauge redundancy as we can still redefine [fi a i, fia2\ T — > 
[fiai,e tXi fia2] T without changing the form of Q a . This 
remaining U(l) gauge redundancy allows us to fix the 
form of Q' a to be Q' a = r'e l ^[cos9,sm9] T , where r' > 0, 
£ [0,2tt) and 9 G [0,tt). 

In addition, we may further demand the mean-field 
Hamiltonian to be time- reversal invariant, which fix <j!> = 
or 7T in the above expression for Q' a . Absorbing the 
sign coming from (f> — n into 9, it thus suffices to take 
Q a = r[l,0] T and Q' a = /[cos 9, sin 9} T , with r, r' > 
and 9 G [0, 2ir). 

In Fig. [18] we plot the results obtained from mini- 
mizing (H t jj') MF with respect to the general ansatze 
parametrized above. From the figure it can be seen that 
states with Q a ^ Q' a emerges only for J > 2 J'. Com- 
paring Fig. [18] with Figs. [T4] and [To] it can be seen that 
the region for which Q a ^ Q' a is deep inside the ferro- 
magnetic phase of the effective spin- 1/2 model and does 
not fit well to the location of the QS phase in the DMRG 
phase diagram Fig. [4] 



However, it should be noted that the low-energy gauge 
structure of an ansatz with 9 =/= is markedly differ- 
ent from that of an ansatz with 9 = 0. Specifically, 
when 9 = 0, the gauge transformation [fi a i,fia2] T ~> 
[fiai,e lXi fia2] T leaves the mean-field ansatz invariant, 
implying that a U(l) gauge field remains gapless in the 
low-energy effective theory. In contrast, when 9^0 there 
is no continuous transformation that leaves the ansatz in- 
variant, hence no gapless gauge field remains in the low- 
energy effective theory. Because of this difference, it is 
conceivable that the energies of mean-field ansatze with 
9 =/= may renormalize differently from those with 9 = 0, 
thus opening the possibility that ansatze with 9^0 may 
become favorable in the region of parameter space that 
corresponds to the QS phase. In what follows, we shall 
stop worrying about the mean-field energetics and in- 
stead focus on whether the low-energy effective theory 
obtained from mean-field ansatze with 9^0 can account 
for the QS phase. 

In the Schwinger fermion decomposition Eq. p]), 



*i±4 



Thus, setting r ^ 



fT r 

r' in the mean-field ansatz will result in a mean-field 
state for which (Si ■ s i+ i^ ^ /Si-s^i^. However, 
from DMRG we know that within numerical accuracy 
^Si-s i+ i^ = ^Sj-Sj_i^ [c.f. Fig. (6)]. Hence, in the 

remaining we shall consider only the case in which r = r' 
but 9^0. 

For illustration, we pick r/t = 1.812, which when 9 = 
corresponds to, e.g., J = 8/3 and J' = 16/3, from our 
previous calculation (a region which from Fig. 



16 



one 

might expect to be proximate to the QS phase), and in- 
troduce an ad hoc value of 9 = 7r/8 to the mean-field 
ansatz. The original mean- field spectrum with 9 = 
and the modified spectrum with 9 = tt/8 are shown in 
Figs. 19(a) and |19(b)[ respectively. Somewhat surpris- 
ingly, the mean-field spectrum is largely unaffected by 
the change of 9. In particular, the flat band at Fermi 
energy continues to appear in the spectrumPS This be- 
havior seems to be a generic feature for this class of 
ansatze. i.e., this flat band exists for general values of 
9 and r/t. However, importantly, the compositions of 
the eigenstates in this flat band is modified, as can be in- 
ferred from the coloring of the band in Fig. |19(b)| (color 
online), which shows that the states near k = ±7r have 
large wavefunction overlaps with the f\ spinons. 



Lifting the degeneracies in the class of extended 
mean-field ansatze 



As before, effects beyond mean-field are expected to 
lift the degeneracy of the flat band. In principle, one can 
apply degenerate perturbation theory as presented in the 
preceding section, but with two significant modifications: 
First, the zero-energy mean-field single-particle state 7jo a 
on site i now has to be constructed from Wannier orbit als. 
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FIG. 19: (color online) Band structures obtained from H MF 
with the additional /-hopping term t' , for r/t = r'/t = 1.812 
and various values of t' and 6. The color encode the eigenvec- 
tor composition of the bands, with red ~ c^, green ~ fik, and 
blue ~ /2fe [the detailed color scheme is shown in (d)] . Fig. (e) 
also defines the convention for the bosonization treatment. 



Since both species of spinons have non-zero Kondo hop- 
pings, 7i0a is no longer locally conserved. Consequently, 
terms of the form 7j 0Q 7j 0ct (i ^ j) can appear in the 
effective Hamiltonian. Second, the constraints enforcing 
terms A and N in H MF now depend crucially on the pre- 
cise manner in which the degeneracy is lifted, and hence 
cannot be left out in the residual interaction H rcs . On 
the technical level, it is challenging to perform the per- 
turbative calculation with the two modifications stated 
above. 

In addition, on the conceptual level, in the QS state 
singularities appear in both the nickel and electron spin- 
spin correlations, the latter of which are absent in the 
other phases in the DMRG phase diagram. Such sin- 
gularities in the electron spin-spin correlation cannot 
be easily captured in the perturbation theory, since the 



Wannier orbitals 7^00 are predominately spinon in char- 
acter. 

Therefore, here we take an alternative approach in 
which an additional ad hoc spinon-spinon hopping term 
t'Yji (fta,afi+i,a,a + h-c) is introduced into the mean- 
field Hamiltonian H MF . Such term can be thought of 
as arising from the spinon mean-field decomposition of 
the nearest neighbor spin-spin interaction term J7i -,-Sj • Sj 
generated by the de gene rate perturbation. 

In Figs. 19(c) and (e) we plot the resulting mean-field 
spectrum for t' = — Q.lt < 0, with r/t = 1.812 and 
9 = 7r/8 as before, from which we see that the resulting 
spectrum now has four Fermi crossings at incommensu- 
rate wavevectors. This four-crossing spectrum appears 
to be a general feature of the mean-field ansatze when 
8^0 and t' ^ 0. i.e., they exists as long as 8 ^ and 
that t' is small but non-zero. However, the Fermi ve- 
locities at the four crossing points will be inverted when 
t! > 0. More importantly, the resulting crossings will 
be predominately spinon in character, while for t' < 
the crossings at k = ±kp2 near ±7r will have a non- 
negligible electron weight (for the parameters used in 
Figs. J19(c)| and [Rl \(ck F2a \lk F2 o a )\ ~ 0.3). We shall 
therefore take t' < and consider the bosonized theory of 
the archetypal band structure shown in Fig. 19(e) Apart 
from the minor complication arising from the matrix ele- 
ments resulting from the compositions of the low-energy 
fermions in terms of the original (c,/i,/2) fermions, the 
bosonized theory of the four-Fermi-crossings band struc- 
ture shown in Fig. 19(e) has been studied extensively by 
Sheng, Motrunich, and Fisher in the context of the so- 
called spin Bose metal (SBM) P We therefore adopt most 
of their notations and keep our account down to the es- 
sentials by referring our readers to Ref. [S]for details. 



C. Bosonization of the mean-field theory 

As in Ref. O we define eight species of low-energy 
fermions ipp aa , one for each Fermi point, in which 
P = R/L = +/— labels the two propagation direc- 
tions, a =t, 4- labels the two spins, and a = 1,2 cor- 
responds the two Fermi wavevectors kp a . The two Fermi 
wavevectors are chosen such that fermions at kp a are 
right-moving, and that |A;f2| > I^fiI (see Fig. 19(e)| for 
illustration). Note that the Fermi wavevectors satisfy the 
relation kFi + &F2 = — 7r/2. 

Next, we bosonize the fermions as follows: 



0i 



oc T] aa e 



(13) 



in which (p aa and 8 aa are bosonic fields that sat- 
isfy [p aa {x),(p b p(x')\ = [8 aa (x),8 b p(x')} = and 
VPa a {x),d b p(x r )] = iTr5 ab 5 al 3<d(x - x') (here ®{x) is the 
Heaviside step function with regularization 9(0) = 1/2), 
while rj aa are the Klein factors satisfying {T] aa ,r] b p} = 
2S ab S a fj. 

As in Ref. |U in addition to the above {1 1, 1 4-, 2 t> 2 1} 
basis for the bosonized fields, it is useful to introduce also 
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the {I/O, la, 2p, 2a} basis and {p+, p—, a+, a—} basis, de- 
fined by the following canonical transformations of the 
6>-ficlds and ^-fields (the transformations for the ^-fields 
are given by replacing every 9 with <p in the equations 
below): 



7af — Pq4. 

hp. - 8<in , s 



(a = 1,2) ; (14) 



(15) 



As customary, we shall refer to fields with index p as the 
charge fields and fields with index a as the spin fields. 

At the level of fermion bilinears, the wavevectors 
±2kp a (a = 1,2), ±7r/2, and ±{kpi — kp2) are brought 
out, and the corresponding bosonized expressions for the 
nickel spin Sk and the electron density 5nk are given by: 



QX 

D 2k Fa 


OC 


e V2i< 


sin(V2ip aa ) , 


(16) 


°2k Fa 


OC 




K " cos{V2ip aa ) , 


(17) 


QZ 

D 2k Fa 


oc 


e V2it 


sm(V29 aa ) , 


(18) 


n 2k Fa 


OC 


e V2U 


>^cos{V29 aa ) ■ 


(19) 



In the absence of any residual interactions, the La- 
grangian density for the bosonized fields is given by: 



act 



v a (d x 



(28) 



It is important to remark that, in our case the bosonized 
fields are the only low-energy degree of freedom remain- 
ing in the theory. More precisely, recall that A and N 
should properly be thought of as fluctuating fictitious 
gauge fields, and that Q a and Q' a should properly be 
thought of as fluctuating bosonic fields. However, for the 
mean-field ansatze that we now consider, all fictitious 
gauge fields have been gapped through the Anderson- 
Higgs mechanism, with various transverse components 
of Q a and Q' a fluctuations serving as the corresponding 
Goldstone boson that are "eaten up." The remaining 
fluctuations of Q a and Q' a are gapped, upon integrat- 
ing out high-energy degrees of freedom if not at the bare 
level. Therefore, contrary to Ref. 6, there is no a priori 
reason for 9 P + to be pinned. 



SI 12 « ( 


. . .)er ie "+e ie '- sin(^ p _ - 


V<y+) 




+ ( 


. . .)f 






(20) 




. . .)e~ i0 "+e i0 "- cos{(p p - - 






+ ( 


. . .)e 


,-ie p+e -ie a - cos (^ p _ _ 


f <P<T+) , 


(21) 


S l/2 « ( 


. . .)e~ ie "+e ie '+ sm((p p _ - 


Per-) 




+ ( 


. . .)f 




~ fa-) , 


(22) 


6n„ /2 oc ( 


. . .)e~ i0 "+e i0 "+ sin(</? p _ - 


fa-) 




+ ( 


. . .)e 




- fa-) ; 


(23) 


QX 

°k F1 -k F2 


oc ( 


..)e l9 "-e- ,9 '+sin(^. 


- - fa+) 






+ (• 


..)e ie "-e ie '+ Bin(vJ p - 


+ fa+) , 


(24) 


k F1 -k F2 


oc ( 


..y 8p -e- t9 °+cos{ip p 


- - fa+) 






+ (■ 


_ye p - e ie„ + cos ( (/?p _ 


+ fa+) , 


(25) 


QZ 

l - f k F i — k F 2 


oc ( 


..y^-e- 10 - sin(^ p 


- - fa-) 






+ (• 


..)e w H-e M -sin(pp_ 


+ fa-) , 


(26) 


0~'H'k F \ — k F 2 


oc ( 


..)e ie "+e- ie "- sin(<p p . 


- ~ fa-) 






+ (• 


..y 9 "+e 10 "- sin(v9 p _ 


+ fa-) ■ 


(27) 



where (. . .) represents various numerical and Klein fac- 
tors, which are not important for our purposes. Also, the 
bosonized expression for electron spin are essentially 
the same as that of the nickel spin except for changes in 
the numerical factors in the (. . .) due to matrix elements. 
As usual, 0_ k = 0\ for Sfc, s^, and 5rik- 



D. Interactions in the bosonized theory 



In the absence of any pinnings of the bosonic 
fields, the low-energy effective theory described by 
Eq. ( p8| is a c = 4 Luttinger liquid. However, as 
the DMRG results show a charge gap in QS phase, 
we shall accept as an empirical matter that 9 p+ is 
pinned, which can happen if the eight-fermion interaction 

i } Rn^RlA^Rn^Ll^Lli^L2^L2i + h - C - « ™s(19 p+ ) 

is sufficiently strong. Moreover, since there are two 
incommensurate Fermi wavevectors kpi and kp2, the 
correlation functions in such theory is expected to ex- 
hibit singularities at multiple incommensurate wavevec- 
tors that are integer combinations of kpi and kp2, as 
well as at wavevectors that arc integer multiples of ir/2 — 
— {kpi — kp 2 ). Since the correlation functions obtained in 
DMRG (Fig. [5| show prominent singularities at merely 
one pair of incommesuate wavevectors k = 7r ± 8, the 



free theory given by Eq. (28) seems to be inconsistent 



with DMRG. Such inconsistency could in-principle be 
accounted for if all the undesired singularities are sup- 
pressed by non-universal amplitudes. However, it is more 
natural to consider scenarios in which some of the bosonic 
fields are pinned by interaction. 

As explained in Ref. [6l assuming that chiral interac- 
tions lead only to velocity renormalizations, the four- 
fermion interactions schematically consists of three parts, 
Ant-4 = W + V p + V a . In terms of the bosonized fields, 
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these read: 
W = cos(2^ p _) 



4w£ 2 (cos(2^ CT _) - f cos(26» CT _ 



wf 2 ( cos(2y) (T _) + f cos(26> ff _) + 2fcos(26» CTH 



(29) 



^ = E^((W 3 -(^«. 

a 

+ ^r((S t 0i P )(0»02„) - (5^i P )(a^ 2p )) , (30) 

^ = E X aa cos{2V2e aa ) + 2A? 2 f cos(20 CT+ ) cos(2<^_) 

a 

A° 



+E#((^-) 2 -(^-) 2 

a 

+ ((^Via)(9 x ^) - (d x la )(d x d 2a )) , (31) 

where w^ 2 an d A^ b (a, & = 1, 2 and a < b; p — p,a) are 
parameters that control the interaction strength. More- 
over, as in Ref. 6 ( the parameters A^ fc satisfy the following 
RG equations at the one-loop level: 



dKa 

dt 



(Kg) 

2ttv„ 



dAf a 



(A 12 ) : 



7r(«i + U 2 ) 



(32) 



such that the only instabilities caused by V a arise from 

K b < o. 

Note that even after incorporating V p , the resulting La- 
grangian density CJ = Co + V p remains quadratic. How- 
ever, scaling dimensions of operators that contain the 
charge fields will be modified. As a result, the terms 
in W will in general acquire scaling dimensions that are 
different from their bare value. Hence, we can consider 
separately the case in which W contains the most rele- 
vant terms and the case in which V a contains the most 
relevant terms. 

If W contains the most relevant terms, then in general 
(p p - is pinned, and correspondingly # p _ is completely 
disordered. From the bosonization formula and from the 
definition of p _ it can be check that the bosonized ex- 
pression of a product of fermions does not contain 8 p - if 
and only if it carries momenta that are multiples ofir/2. 
Consequently, the pinning of </? p _ will kill the singulari- 
ties of any correlation functions located at incommensu- 
rate wavevectors. This is inconsistent with the DMRG 
results. Hence, we conclude that in the bosonized theory 
of the QS phase W must be irrelevant. 

Next we consider the case in which V a contains 
marginally relevant terms. In the simplest scenarios, the 
terms with which A^ fc < contain mutually commuting 
variables. There are four such cases: 

(i) Afj > 0, A£ 2 > 0, \<[ 2 < ; 

(ii) < 0, A§2 < 0, \\ 2 > ; 

(iii) A?! < 0, Af 2 > 0, A? 2 > ; 



(iv) AJj > 0, Afa < 0, \1 2 > . 

In cases ^ and ([TTJ> two mutually commuting spin fields 
are pinned. Consequently, all spin-spin correlations must 
be void of singularities. Thus, these cases are inconsistent 
with the DMRG results. In cases ( [iii] ) and ( |rvj ), one spin 
field associated with wavevector kp a is pinned. Thus, 
at the level of fermion bilinears, the only singularities 
in the spin-spin correlations that survive are the ones at 
±2k Fa ,, where a' ^ a [c.f. Eqs. @-([27|]. These singu- 
larities can be identified with the singularities present in 
DMRG QS phase. Moreover, recall that in the DMRG 
results there are broad peaks in the nickel (but not elec- 
tron) spin-spin correlation located roughly at wavevec- 
tors 7r± /c s i n g, where ±/c s j ng are the wavevectors of which 
the prominent singularities are seen. These broad peaks 
can be interpreted as the remnants of the singularities 
at ±2kp a after the associated spin field is pinned. Re- 
call that the Fermi points at ±fc_p 2 have non-negligible 
electron character while the Fermi points at ikFi are 
predominately spinon in character. Since in the DMRG 
the singularities are observed in both the nickel and the 
electron spin-spin correlation, while the broad peaks are 
observed only in the nickel spin-spin correlation, we may 
identify the state obtained in DMRG with scenario 
Since we also assume that # p + is pinned, the resulting 
state would be a c = 2 Luttinger liquid with one charge 
mode and one spin mode ("C1S1"). 

One potential objection to this identification is that 
it implies that the electron density-density correlation 
function (Sn^^du^) possesses singularities at ±2kpi and 
±2kp 2l which are not observed. However, it is known 
that when the charge fluctuation is reintroduced to the 
SBM, as the system approaches the Mott transition and 
when it is in the insulating phase, the non-universal am- 
plitudes in the density-density correlation can be suffi- 
ciently small that in numerics it can appear smoothP^ 
The situation that we found in the DMRG study of the 
present model may correspond to such a situation. 

More generally, one might pose the question of whether 
the DMRG results can be explained by other combina- 
tions of pinned fields, which can arise from higher-order 
interactions. Here we briefly consider such possibilities. 

Since we assume that 9 p+ is pinned, the {p+, p—} 
basis is the appropriate basis to describe the charge 
sector. The only question left for the charge sector is 
whether ip p _ or # p _ can be pinned. As already men- 
tioned, the pinning of y? p _ would kill all correlations at 
incommensurate wavevectors and hence is inconsistent 
with the DMRG results. As for 6 p _, one can check that 
it carries an incommensurate momentum kp\ — kp 2 (i.e., 
# p _ — > (9 p _ + kpi — kp 2 under the translation i i + 1), 
and hence cannot be pinned. Thus, unless one appeals to 
vanishingly small non-universal amplitudes, exactly one 
spin degree of freedom must be gapped to produce the 
single pair of prominent singularities in the spin-spin cor- 
relation observed in DMRG. Therefore, we are left with 
a "C1S1" case similar to the one we analyzed, except 
that the pinned spin field may more generally be a linear 
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combination of the la and 2a fields. 



VI. DISCUSSIONS AND CONCLUSIONS 

In this manuscript we considered a Id underscreened 
Kondo chain with alternating spin-one and electron sites, 
which in addition to the familiar electron hopping and 
Kondo term also contained a spin-dependent hopping 
term. We analyzed the model numerically using DMRG 
and found that the phase diagram consists of a ferro- 
magnetic (F) phase, a quasi-long-range antiferromagnetic 
(QAF) phase, a charge-density (CD) ordered phase, and, 
importantly, a quasi-long-range spiral (QS) phase, in 
which singularities in the spin-spin correlation are ob- 
served at incommensurate wavevectors. 

To interpret the DMRG results, we introduced a slave- 
particle representation of the spin-one spins, from which 
we obtained a mean-field Hamiltonian. Taking the sim- 
plest class of mean-field ansatze and making use of 
second-order degenerate perturbation theory, we were 
able to map the mean-field Hamiltonian to an effective 
spin-1/2 model, from which follows a quantum mean- 
field phase diagram that resembles the one obtained from 
DMRG, with the exception of the QS phase. We then fo- 
cused on the QS phase, by extending our class of mean- 
field ansatze and considering the bosonized interacting 
theory arose from them. By considering various pos- 
sible interactions at the four-fermion level, we argued 
that the QS phase is best reproduced by an interacting 
bosonized theory with one charge and one spin degree 
of freedom ("C1S1"), in which the spin fields carry an 
incommensurate wavevector. This particular interacting 
bosonized theory is proximate to the spin Bose metal 
phase proposed by Sheng, Motrunich, and Fished. Con- 
sequently, our results point to a possible route to obtain 
gapless quantum spin liquids with singularities at incom- 
mensurate wave vecto rs that does not involve the "ring- 
exchange" terms f^l and may pave the way for realizing 
such states in Id as well as constructing similar states in 
higher dimensions. 

While we believe that the current work presents a co- 
herent and comprehensive analysis of the system we con- 
sidered, we also mention a few opportunities for further 
study. On the numerical side, it is clearly desirable to ob- 
tain the central charge of the QS phase from DMRG and 
check against the c = 2 prediction from mean-field the- 
ory. This requires either further computation resources 
or improvements of the methodology. Moreover, it would 
be helpful to check the projected energetics of the mean- 
field ansatze we proposed, and compare them against the 
energy obtained from DMRG, via variation Monte Carlo 



(VMC). Unfortunately, because the slave-particle repre- 
sentation we used requires additional constraints on top 
of the local conservation of spinon number (c.f. Eq. |3j, 
such a calculation is challenging. On the analytic side, it 
would be desirable to develop a method that can directly 
obtain the CD phase from the slave-particle representa- 
tion without appealing to an indirectly argument based 
on symmetry. Moreover, we note that there is an alter- 
native slave-pa rticle representation of the spin-one spins 
in the literatur e ! 25 1 26 1 , and it would be good to check if re- 
sults similar to those we obtained can be derived in this 
alternative slave-particle representation. 

Finally, we return to the nickel valence controversy 
that motivated this work, which in the simplest case can 
be considered as a three-dimensional (3d) extension of 
the present model. From our mean-field analysis and 
from general arguments on Kondo problems we are quite 
convinced that the low-energy physics of such 3d model 
will still be described by an effective spin-1/2 model, 
which tends to exhibit static spin orders. However, the 
dimer phase in the effective spin-1/2 model, which we ar- 
gue to correspond to the oxygen-centered charge-density 
order, is a special feature in Id, and is unlikely to per- 
sists to 3d. In any case, our DMRG phase diagram shows 
no sign of nickel-centered charge order (or its symmetry- 
equivalent order — we do not insist upon any significant 
charge accumulation) , which has been interpreted as due 
to local Kondo singlet formation. Thus, the direct exten- 
sion of our model to 3d is unlikely to exhibit any charge 
order, whether it is oxygen-centered or nickel-centered. 
In our opinion this points to the necessity of including 
explicit coupling to oxygen motions in the mechanism of 
charge ordering. It would be interesting to explore this 
directly by explicitly including lattice degrees of freedom 
in a model similar to ours. We leave this topic for future 
study. 
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